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EXECUTIVE  SUMMARY 


V 

We  have  applied  a  system  modelling  approach  to  the  problem 
of  projecting  the  performance  of  an  enhanced  NAVSPASUR  system. 
We  have  modelled  the  performance  characteristics  of  the  current 
system  as  well  as  the  out-of-plane  station  (OOPS)  and  enhanced 
signal  processing  options  described  in  the  NAVSPASUR  System 
Development  Plan,  Phase  II.  The  objective  of  this  modelling 
process  was  to  evaluate  the  ability  of  NAVSPASUR  to  provide 
satellite  position  and  velocity  measurements  to  the  accuracy 
required  for  first  pass  orbit  determination  (FPOD) .  The  FPOD 
position  and  velocity  requirements  have  been  taken  to  be  1.67  km 
and  2  km/min,  respectively,  as  specified  in  the  NAVSPASUR  System 
Development  Plan,  Phase  I.  |?  >4  /  lutyi- 1  ^  ^ 

The  major  results  of  this  effort,  described  in  detail  in 
this  report,  are  summarized  below. 

1.  The  current  system  appears  to  be  performing  near 
optimal  levels  in  terms  of  the  resultant  position  and 
velocity  uncertainties  for  single  pass  data.  For 
satellites  with  a  1  meter2  effective  radar  cross- 
section,  accurate  position  determinations  are  limited 
to  objects  crossing  the  fence  over  the  continental 
U.S.  (CONUS)  at  altitudes  below  about  3000  km. 
Velocity  determinations  are  insufficient  to  meet  FPOD 
requirements  at  virtually  any  fence  crossing 
location.  Typical  velocity  errors  range  between  3  and 
18  km/minute  over  much  of  the  fence  at  altitudes 
between  200  and  2000  km.  There  is  a  conspicuous 
degradation  in  low  altitude  system  performance  to  the 
west  of  longitude  100°  W,  with  velocity  errors  up  to 
35  km/min  and  positional  errors  in  excess  of  3  km  at 
altitudes  below  200  km.  For  altitudes  greater  than 
2000  km,  velocity  errors  generally  exceed  20  km/min 
across  the  full  extent  of  the  fence. 

2.  The  improvement  of  Doppler  tracking  capability  to  *1 
Hz  and  the  addition  of  Doppler  rate  (chirp)  tracking 
to  ±4  Hz/sec  results  in  typical  improvements  in 
velocity  determination  of  30  -  50%  over  most  of  the 
fence,  while  positional  accuracy  remains  effectively 
unchanged.  The  effect  of  the  western  "hole”  is 
somewhat  diminished,  but  not  entirely  eliminated. 
However,  even  at  this  level  of  performance,  the 
system  is  incapable  of  meeting  the  specified  velocity 
requirements  at  most  fence  crossing  locations. 

The  addition  of  three  new  out-of-plane  receiving 
stations  (OOPS)  with  phase  and  phase  rate  measurement 
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capability  to  the  current  array  (assuming  ±10  Hz 
Doppler  processing  and  no  chirp  measurement)  produces 
two  significant  improvements  over  the  present  system. 
First,  there  is  an  overall  improvement  of  more  than 
an  order  of  magnitude  in  both  position  and  velocity 
accuracy  for  all  illuminated  objects.  Secondly,  with 
a  judicious  choice  of  the  OOPS  locations,  the  western 
"hole"  in  the  fence  coverage  is  virtually  eliminated. 
For  this  system,  velocity  errors  are  less  than  2 
km/min  for  virtually  all  CONUS  fence  crossings  at 
altitudes  up  to  3000  km.  Position  errors  meet  the 
stated  requirement  for  all  orbits  investigated,  with 
the  exception  of  100  km  orbits  in  the  far  west. 

4.  If  ±10  Hz/sec  chirp  processing  is  added  to  the  above 
system,  only  marginal  improvement  is  realized. 
However,  if  the  Doppler  and  chirp  measurement 
capabilities  are  enhanced  to  achieve  accuracies  of  ±1 
Hz  and  ±4  Hz/sec,  respectively,  an  additional  factor 
of  3  improvement  is  attained.  For  this  case, 
velocity  errors  are  typically  between  0.1  and  0.6 
km/min  up  to  altitudes  of  4000  km  for  all  CONUS  fence 
crossings,  except  for  100  km  altitudes  at  the  edges 
of  CONUS.  Performance  up  to  3°  longitude  off  the 
coast  is  within  1  km/min  at  altitudes  of  500-4000  km. 
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The  addition  of  3  OOPS  which  provide  only  Doppler 
and  chirp  measurements  could  meet  the  velocity 
requirements  for  almost  all  CONUS  fence  crossings  at 
altitudes  up  to  1600  km,  assuming  Doppler  and  chirp 
measurement  accuracies  of  ±10  Hz  and  ±10  Hz/sec, 
respectively.  If  measurement  accuracies  of  ±1  Hz  and 
±4  Hz/sec  can  be  achieved,  velocity  errors  are 
reduced  to  less  than  1  km/min  up  to  4000  km 
altitudes.  For  both  cases  the  position  errors  are 
within  the  required  accuracy  for  altitudes  up  to  4000 
km,  with  the  exception  of  100  km  orbits  in  the  west. 

The  results  of  our  analysis  indicate  that  the  target 
position  and  velocity  requirements  can  be  met  with  an 


enhanced  NAVSPASUR  system  which  appears  to  be  both ^ — 
technically  and  economically  feasible.  However,  we  n  or_ 
have  not  at  this  point  addressed  the  specifics  of  the  A41 
design,  software  and  hardware  requirements,  data' 
transfer  protocols,  or  data  analysis  algorithms  which 
would  be  required  to  implement  these  enhancements.  ^ 
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I.  Introduction 

We  have  developed  a  set  of  algorithms  which  are  designed  to 
accurately  model  the  current  performance  characteristics  of  the 
NAVSPASUR  full -Doppler  system  both  in  terms  of  amplitude  response 
and  inherent  phase  accuracy.  We  have  applied  these  algorithms  to 
a  family  of  satellite  orbits  to  generate  model  NAVSPASUR  phase 
and  amplitude  data  sets.  These  model  data  were  then  input  into  a 
nonlinear  least  squares  fitting  routine  used  to  estimate  the 
likely  position  and  velocity  uncertainties  as  a  function  of 
satellite  radar  cross-section,  altitude,  longitude  at  fence 
crossing,  and  orbital  inclination.  In  addition  to  performing 
this  analysis  for  the  current  system  configuration,  we  have 
extrapolated  current  receiver  site  performance  to  proposed 
additional  out-of-plane  stations  (OOPS) ,  and  have  evaluated  the 
potential  impact  such  stations  would  have  on  overall  system 
performance. 

We  have  also  examined  the  effect  of  improved  Doppler  signal 
processing  for  both  the  current  system  and  an  expanded  system 
including  three  OOPS.  The  impact  of  adding  chirp  processing  to 
both  the  current  and  expanded  systems  was  examined  as  well. 
Finally,  we  explored  the  consequences  of  adding  OOPS  which 
provide  Doppler  and  chirp  tracking  only,  with  no  phase 
measurement  capability. 

We  describe  below  the  basic  assumptions  we  have  adopted  for 
this  investigation.  The  details  of  the  data  modelling  and  error 
analysis  methods  are  described  in  Sections  II  and  ill,  while  the 
analysis  of  our  results  is  contained  in  Section  IV.  In  Section  V 
we  discuss  the  areas  in  which  further  research  is  either  ongoing 
or  required.  Finally,  in  Section  VI  we  summarize  the  results  of 
the  current  research  effort  and  discuss  their  implications  for 
future  changes  to  the  NAVSPASUR  system. 

a).  Approach  and  Basic  Assumptions 

The  primary  objective  of  this  project  was  to  predict  the 
system  performance  of  an  enhanced  NAVSPASUR  system  as  envisioned 
in  the  options  presented  in  the  System  Development  Plan,  Phase 
II.  For  a  system  as  complex  as  a  multi-station  bi-static  radar, 
the  net  uncertainties  in  the  satellite  position  and  velocity 
which  result  from  known  uncertainties  in  the  measured  quantity, 
antenna  phase,  depend  nonlinearly  on  both  the  phase  measurement 
errors  and  the  transmitter-satellite-receiver  geometric 
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relationships.  For  such  a  system  the  most  straightforward  method 
of  predicting  system  performance  is  a  modelling  approach. 

For  the  current  system  we  have  operational  data  available 
against  which  model  data  can  be  compared.  Given  the  availability 
of  this  data  one  might  expect  that  the  analysis  of  current  system 
performance  should  proceed  solely  on  the  basis  of  real  data. 
While  we  concede  the  value  in  such  an  approach,  we  have  opted  to 
base  our  analysis  on  a  system  model  approach  for  several  reasons. 

First,  it  is  clear  that  the  real  data  contain  systematic 
errors  which  are  at  present  poorly  understood  and  which  may  be  in 
part  correctable.  These  include  both  phase  calibration  errors 
and  clock  offset  errors,  among  others.  Their  inclusion  in  the 
current  system  analysis  (by  virtue  of  using  real  data)  but  not  in 
the  enhanced  system  analysis  would  make  a  comparision  between  the 
different  systems  difficult.  Secondly,  the  necessity  to  include 
clock  offsets  in  a  least  squares  solution  using  real  data  would 
double  the  number  of  parameters  to  be  solved  for  and  would 
increase  the  computational  load  significantly.  Finally,  the 
modelling  approach  allows  us  to  specify  the  target  size  and 
satellite  orbital  parameters,  so  that  we  can  investigate  the 
system  performance  over  the  full  range  of  satellite  orbits  in 
which  we  are  interested. 

Although  our  formal  analysis  of  the  current  system 
performance  is  limited  to  modelled  data,  we  have  compared  our 
results  to  the  daily  aggregate  system  performance  reports 
compiled  at  Dahlgren.  We  find  that  the  results  we  obtain  from 
our  system  model  are  consistent  with  these  summary  reports, 
indicating  that  the  approach  we  have  adopted  is  basically  sound. 

In  our  analysis  we  have  made  several  assumptions  which  we 
enumerate  here,  deferring  a  discussion  of  their  significance 
until  a  later  point. 

1.)  The  current  system  data  acquisition  rate  (nominally 
1/55  second)  and  phase  data  quantization  factor  (1/64 
rotation,  or  5? 6)  would  be  retained  in  any  future 
operations  as  studied  in  this  report. 

2)  .  The  operation  of  any  additional  OOPS  will  be  such 

that  the  phase  difference  measurement  errors  would  be 
qualitatively  the  same  and  of  the  same  RMS  magnitude 
as  those  for  the  current  in-plane  stations. 

3)  .  The  measurement  errors  inherent  in  the  determination 

of  the  individual  antenna  array  phases  are  randomly 
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distributed  and  vary  smoothly  with  received  signal 
amplitude  according  to  the  error  model  developed  in 
the  report  The  Determination  of  NAVSPASUR  Phase  Errors 
(hereafter  Paper  I) .  Any  systematic  errors  present 
in  the  system  are  assumed  to  be  either  correctable  or 
of  sufficiently  small  magnitude  so  that  their 
presence  does  not  significantly  alter  the  results  we 
obtain. 

4) .  A  Doppler  search  and  tracking  accuracy  of  ±1  Hz  is 
attainable  with  available  hardware  using  standard 
signal  processing  techniques. 

II.  Data  Modelling  Technique 

The  primary  objective  in  generating  the  model  data  was  to 
produce  data  scans  which  closely  resembled,  in  qualitative  and 
quantitative  terms,  the  real  data  acquired  by  the  NAVSPASUR 
system  for  satellites  of  similar  cross-section  and  orbits.  The 
steps  involved  in  this  process  were: 

1)  .  Generate  a  family  of  satellite  orbital  elements 

corresponding  to  the  cases  to  be  investigated. 

2)  .  Develop  an  analytic  model  for  transmitter  and 

receiver  beam  patterns  which  produce  an  acceptable 
received  amplitude  response  as  a  function  of 
satellite  position. 

3) .  For  each  data  collection  time  interval,  generate  the 

ideal  phase  differences  expected  for  each  independent 
baseline  pair  at  each  receiving  station. 

4) .  Introduce  random  noise  on  the  ideal  phase  differences 

using  the  error  model  obtained  from  the  phase 
difference  error  analysis  of  Paper  I. 

5)  .  Calculate  the  ideal  Doppler  frequency  and  chirp  for 

each  receiver-transmitter  combination. 

6)  .  Introduce  random  errors  on  the  ideal  Doppler  and 

chirp  values  according  to  the  assumed  measurement 
accuracy . 

The  details  of  the  data  modelling  algorithm  are  outlined  in 
the  following  sections. 
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a).  Satellite  orbital  elements 


We  have  investigated  the  system  response  to  satellite  passes 
at  altitudes  of  100,  200,  300,  500,  700,  1000,  1300,  1600,  2000, 
2500,  3000,  3500,  and  4000  kilometers  which  penetrate  the  fence 
at  longitudes  ranging  from  7  5°  to  12  0°  west  longitude  in  5° 
increments.  The  full  altitude-longitude  grid  was  explored  for 
satellites  with  an  assumed  orbital  inclination  of  85°.  This 
inclination  represents  a  relatively  unfavorable  case  in  that 
these  satellites  are  moving  nearly  perpendicular  to  the  NAVSPASUR 
great  circle  and  therefore  have  a  nearly  minimum  beam  crossing 
time.  In  order  to  check  the  sensitivity  of  our  results  to  the 
assumed  orbital  inclination,  we  investigated  a  subset  of  the 
altitude-longitude  grid  at  orbital  inclinations  of  40°  and  60°. 
As  expected,  the  results  obtained  at  these  smaller  orbital 
inclinations  were  marginally  better  than  those  at  an  85° 
inclination. 

All  orbits  were  assumed  to  be  circular  for  computational 
convenience.  This  restriction  is  not  expected  to  have  any 
significant  effect  on  our  results,  at  least  for  objects  whose 
orbital  eccentricities  are  not  excessive.  However,  in  the  case 
of  satellites  with  highly  elliptic  orbits  crossing  the  fence  near 
perigee,  some  degradation  of  system  performance  from  the  levels 
predicted  herein  should  be  expected,  given  their  relatively 
higher  orbital  velocities  and  correspondingly  shorter  fence 
crossing  times.  Future  extensions  of  this  work  will  address  the 
case  of  highly  elliptic  orbits,  since  some  objects  of  interest  do 
in  fact  have  such  orbits. 

For  each  orbit  investigated,  we  calculate  the  position  and 
velocity  of  the  satellite  at  the  epoch  of  fence  crossing  and 
express  these  in  a  rotating,  geocentric  coordinate  system  as  used 
by  NAVSPASUR.  These  positions  and  velocities,  together  with  the 
assumed  radar  cross-section,  provide  the  basic  input  to  the  data 
generation  algorithms. 


b).  Transmitter  and  receiver  beam  patterns 

The  transmitter  and  receiver  gains  as  a  function  of 
satellite  position  were  determined  from  analytic  formulae  which 
express  the  respective  antenna  gains  as  a  function  of  north-south 
and  east-west  direction  angles.  The  analytic  formulae  used 
varied  depending  on  the  site. 


For  each  of  the  transmitter  sites,  the  E-W  beam  patterns 
were  determined  from  an  9th  order  polynomial  fit  to  the  beam 
pattern  given  in  the  report  Radiation  Pattern  Calculation  of  NAVSPASUR 
Transmitter  Element ,  by  Dr.  Steven  L.  Berg  of  Interferometries,  Inc. 
This  beam  pattern  was  obtained  from  calculations  for  an  inverted- 
V  (arrowhead)  dipole  antenna  mounted  horizontally  above  a  finite 
ground  screen  of  the  proper  dimensions.  This  pattern  effectively 
meets  the  E-W  design  specifications  established  for  the 
performance  of  the  transmitter  antenna  arrays  and  is  believed  to 
accurately  represent  the  current  performance  of  the  transmitter 
sites. 

The  N-S  antenna  pattern  for  the  Kickapoo  transmitter  was 
modelled  as  a  Gaussian  with  a  full  width  to  half  power  in  the  far 
field  of  0  °  042 ,  which  is  the  design  specification  beam  width. 
Since  the  far  field  for  the  Kickapoo  antenna  array  begins  at  a 
distance  of  approximately  15,000  km  we  have  included  a  term  in 
the  N-S  beam  pattern  to  compensate  for  beam  broadening  due  to  the 
target  being  in  the  near  field  of  the  antenna  array.  We  have 
compensated  for  this  beam  broadening  by  introducing  an  additive 
term  to  the  Gaussian  full  width  to  half  maximum  of  farray/3<fsat  where 
/array  is  the  overall  N-S  length  of  the  Kickapoo  transmitter  array 
and  efsal  is  the  transmitter-satellite  distance.  The  factor  of 
three  which  appears  in  the  correction  term  was  determined 
empirically  by  comparing  the  model  data  to  real  NAVSPASUR  data. 

The  N-S  patterns  of  the  two  remaining  transmitters  were 
modelled  as  slot  antennas  of  the  appropriate  N-S  length.  For  all 
transmitters  the  gains  were  normalized  so  that  the  integrated 
power  over  the  upper  half -plane  was  equal  to  the  total  radiated 
transmitter  power. 

The  E-W  antenna  patterns  for  the  receiver  sites  were 
modelled  using  the  analytic  formula  for  a  horizontal  dipole 
antenna  0.28  wavelengths  above  an  infinite  ground  plane  of 
infinite  conductivity.  For  the  real  receiver  arrays,  the  dipoles 
are  located  0.322  wavelengths  above  the  ground  screen.  The 
somewhat  lower  height  we  use  was  chosen  to  compensate  for  the 
beam  narrowing  due  to  the  finite  ground  screen.  Using  our 
adopted  model  we  obtain  an  E-W  full  width  , to  half  power  which  is 
the  same  as  the  design  specifications  for  the  receiver  arrays. 

The  receiver  N-S  patterns  were  modelled  as  slot  antennas  of 
the  appropriate  length.  The  resultant  N-S  patterns  agree  well 
with  the  design  specifications  for  both  the  high-altitude  and 
low-altitude  stations.  The  receiver  patterns  were  normalized  to 
produce  unity  gain  integrated  over  the  upper  half-plane. 
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c).  Generation  of  the  ideal  amplitude  and  phase  data 


For  each  of  the  test  orbits  described  in  Section  Ila,  we 
generated  a  set  of  model  data  scans.  A  separate  data  scan  was 
generated  for  each  receiver/transmitter  combination  for  which  the 
satellite  was  at  least  2°  above  the  horizon  at  both  the  receiver 
and  transmitter,  and  for  which  the  received  signal  amplitude  was 
-152  dBm  or  greater  for  at  least  one  data  collection  interval 
(nominally  1/55  sec).  The  amplitude  cutoff  of  -152  dBm 
represents  a  signal-to-noise  ratio  of  approximately  1:1  in  the 
full-Doppler  mode,  which  is  a  reasonable  detection  threshold 
based  on  an  inspection  of  full-Doppler  mode  NAVSPASUR  data. 

For  each  combination  of  transmitter/receiver  pairs  which  met 
the  elevation  test,  we  began  calculating  the  expected  received 
amplitude  40  time  intervals  (nominally  40/55  sec)  prior  to  fence 
crossing.  The  amplitudes  were  calculated  based  on  the  satellite 
position  at  each  time  interval,  the  assumed  transmitter  and 
receiver  beam  patterns,  and  the  assumed  effective  radar  cross 
section.  If  no  amplitudes  greater  than  the  cutoff  were  found 
after  110  time  interval  increments,  the  calculation  was 
terminated  and  no  data  scan  was  written. 

If  an  amplitude  greater  than  the  cutoff  is  found,  we  write  a 
data  scan  header,  compute  and  record  the  received  Doppler 
frequency  (to  the  nearest  1  Hz)  and  the  data  start  time,  and 
begin  calculating  the  expected  ideal  phases.  A  separate  data 
line  is  calculated  and  written  for  each  time  interval  during 
which  the  amplitude  remains  above  the  cutoff.  Once  the  received 
amplitude  falls  below  the  cutoff  or  the  number  of  data  lines 
written  equals  55,  the  data  scan  is  terminated  and  we  proceed  to 
the  next  transmitter/ receiver  pair.  No  secondary  scans  for  the 
same  receiver/transmitter  pair  are  ever  generated,  even  if  the 
signal  remains  above  the  cutoff  after  55  data  lines  are  written. 

For  each  data  line  whose  amplitude  exceeds  the  cutoff,  we 
calculate  the  expected  antenna  phase  difference,  taken  with 
respect  to  the  designated  reference  antenna,  for  each  antenna 
array  at  the  receiving  site.  This  difference  is  simply  the 
difference  in  path  length  from  the  satellite  to  the  respective 
antenna  arrays,  expressed  in  wavelengths.  For  a  receiver  site 
with  n  antenna  arrays,  we  generate  n-1  phase  differences. 
Implicit  in  this  method  is  the  reasbnable  assumption  that 
atmospheric  and  ionospheric  effects  introduce  no  additional 
differential  path  length  changes. 

Each  data  line  written  to  the  data  scan  consists  of  a 
received  amplitude  for  that  time  interval  and  n-1  ideal  phase 
differences.  Each  data  scan  consists  of  between  1  and  55  data 
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lines  for  a  given  transmitter/receiver  combination,  each  with  a 
received  amplitude  greater  than  or  equal  to  -152  dBm. 

d).  The  addition  of  errors  to  the  data 

Once  an  ideal  data  scan  has  been  completed,  we  introduce 
errors  onto  the  ideal  phase  differences  to  simulate  the  errors 
expected  during  normal  system  operation.  The  phase  difference 
error  model  we  have  adopted  for  this  work  is  the  error  model 
produced  by  Interferometries,  Inc.,  which  is  described  in  Paper 
I.  Using  this  model  we  obtain  errors  which  are  randomly 
distributed,  are  zero  in  the  mean,  and  whose  RMS  magnitude  is  a 
function  of  the  received  signal  amplitude. 

For  each  phase  difference  datum  we  use  the  adopted  error 
model  to  generate  a  random  error  based  on  the  received  amplitude 
associated  with  that  datum  and  add  this  error  to  the  ideal  phase. 
This  process  is  repeated  for  all  the  data  points  in  the  scan. 

We  do  not  introduce  errors  to  the  calculated  Doppler 
frequency  and  chirp  at  the  time  the  data  scans  are  written.  We 
defer  adding  errors  to  the  Doppler  data  until  the  least  squares 
fitting  step,  allowing  us  to  use  the  same  input  data  set  for  all 
Doppler  accuracy  cases.  When  the  Doppler  errors  are  added,  they 
are  random  errors  with  zero  mean,  as  in  the  case  of  the  phase 
errors  previously  discussed. 

e).  OOPS  location  and  performance  specifications 

A  major  objective  of  this  work  was  to  estimate  the 
performance  improvements  that  could  be  realized  by  the  addition 
of  receiver  stations  located  out  of  the  NAVSPASUR  great  circle 
plane.  In  the  current  work  we  have  considered  an  enhanced  system 
which  includes  three  OOPS.  We  believe  this  to  be  the  minimum 
number  of  additional  stations  necessary  to  meet  the  target 
position  and  velocity  accuracies  from  coast  to  coast.  We  base 
this  conclusion  on  a  preliminary  aperture  systhesis  analysis 
using  standard  radio  astronomical  techniques,  from  which  we 
obtained  the  nominal  positions  for  the  three  OOPS  considered  in 
this  study.  The  geographic  locations  of  the  three  OOPS  are 
depicted  on  the  map  given  in  Figure  1,  and  are  as  follows: 

N.  Latitude  W.  Longitude 

Northeast  OOPS  35.09  92.59 

Central  OOPS  42.55  97.54 

Northwest  OOPS  37.17  103.38 
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The  results  described  in  this  report  detail  the  type  of 
performance  enchancement  that  can  be  obtained  with  the  addition 
of  multiple  OOPS  stations.  However,  we  do  not  believe,  nor  do  we 
represent,  that  the  configuration  used  herein  is  the  only,  or 
even  the  optimal,  configuration  for  an  OOPS  enhancement.  We  are 
continuing  to  study  the  range  of  possible  station  locations,  as 
well  as  the  operational  performance  requirements  for  the  OOPS. 
The  results  of  these  further  studies  will  be  provided  when 
available. 

The  results  we  derive  in  this  paper  are  based  on  the 
assumption  that  each  OOPS  antenna  array  could  produce  an 
effective  antenna  gain  of  35  dB  in  the  direction  of  the 
satellite.  Preliminary  investigations  indicate  that  similar 
accuracies  are  obtainable  with  lower-gain  OOPS,  with  the  tradeoff 
that  FPOD  capability  would  be  limited  to  lower-altitude  orbits. 
Since  the  question  of  required  gain  is  intimately  related  to  the 
number  and  location  of  the  OOPS,  we  have  deferred  a  detailed 
analysis  of  the  gain  requirements  until  a  later  report,  at  which 
time  we  will  address  the  conceptual  design  of  the  OOPS. 

For  the  models  which  assume  OOPS  with  phase  measurement 
capability,  the  OOPS  stations  consist  of  five  antenna  arrays 
arranged  in  a  filled-cross  configuration  with  baselines  aligned 
in  the  N-S  and  E-W  directions.  The  spacing  between  antenna 
arrays  in  both  directions  is  taken  to  be  5  km. 

III.  Least  Squares  Analysis  of  System  Performance 

The  expected  errors  in  the  derived  position  and  velocity  for 
a  given  satellite  orbit  are  estimated  using  a  recursive  nonlinear 
least  squares  fitting  and  error  estimation  routine  with  data 
weighting,  similar  to  that  described  in  Paper  I.  While  the  least 
squares  technique  is  not  a  viable  algorithm  for  routine 
processing  of  real  NAVSPASUR  data  (due  to  the  requirement  that 
one  have  a  good  a  priori  estimate  of  the  satellite  position  and 
velocity),  in  the  case  of  model  data,  where  the  a  priori  position 
is  known,  it  does  provide  a  valid  means  to  estimate  the  accuracy 
which  can  be  achieved  by  the  system.  For  the  present 
application,  the  least  squares  fitting  is  applied  to  all  data 
scans  for  a  particular  satellite  pass  simultaneously  -  i.e. ,  all 
available  combinations  of  receiving  and  transmitting  sites  are 
used  to  obtain  estimates  for  the  satellite  position  and  velocity 
and  the  likely  error  in  these  estimates. 

The  least  squares  routine  yields  estimates  for  six 
parameters  -  three  position  and  three  velocity  coordinates  at  the 
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epoch  of  fence  crossing,  from  which  the  satellite  orbital 
elements  can  be  obtained.  Since  we  deal  here  with  modelled  as 
opposed  to  real  data,  we  do  not  need  to  solve  for  any  additional 
parameters  such  as  propagation  effects  or  relative  errors  in  the 
clock  times  between  the  different  receiving  sites.  In  fact,  the 
data  generation  and  least  squares  fitting  programs  are  set  up  to 
do  all  calculations  in  "satellite"  time,  so  that  light  travel 
times  are  explicitly  corrected  for  in  the  data  modelling  and 
fitting  routines  by  treating  them  identically  in  both  cases. 

The  coordinate  system  which  we  adopt  in  the  multi-station 
least  squares  fitting  process  is  based  on  the  NAVSPASUR  great 
circle,  so  that  the  x-y  coordinate  plane  is  coincident  with  the 
NAVSPASUR  great  circle  plane  and  the  z  axis  is  normal  to  the 
great  circle  plane.  Thus  the  coordinate  system  is  geocentric  and 
rotates  with  the  rotating  earth.  This  choice  of  coordinate 
system  allows  us  to  determine  the  position  and  velocity  errors  in 
a  coordinate  system  which  is  related  to  the  special  geometry  of 
the  NAVSPASUR  system.  In  this  coordinate  system  the  relative 
magnitudes  of  the  in-plane  and  plane-normal  errors  provide 
insight  into  the  limitations  imposed  by  the  current  geometry  as 
well  as  indicating  the  favored  geometry  for  any  system 
improvements.  The  diagram  shown  in  Figure  2  shows  the  relative 
orientation  of  the  NAVSPASUR  great  circle  coordinate  system  in 
the  rotating  geocentric  reference  frame. 

In  addition  to  generating  phase  data,  we  also  calculate  the 
expected  Doppler  and  chirp  for  each  transmitter-receiver 
combination.  We  have  modified  the  least  squares  program  to  make 
use  of  this  information  in  obtaining  a  best  estimate  of  position 
and  velocity.  When  Doppler  and  chirp  data  are  included  in  a 
given  test  case,  we  introduce  random  errors  on  their  nominal 
values  according  to  their  assumed  accuracy.  As  with  the  phase 
data,  these  data  are  appropriately  weighted  by  their  RMS  error  in 
the  least  squares  fitting  routine. 

Once  the  least  squares  solutions  to  the  satellite  position 
and  velocity  have  been  obtained,  we  also  have  available  the  error 
matrix  F~l  (see  Paper  I)  .  The  diagonal  terms  of  the  error  matrix 
are  the  variances  for  each  of  the  parameters,  while  the  off- 
diagonal  elements  are  the  covariance  terms,  which  give  the 
correlations  between  the  errors  in  the  parameters.  From  the 
error  matrix  we  can  calculate  the  expected  error  in  each 
component  of  the  position  and  velocity,  as  well  as  the  expected 
total  (3-dimensional)  error  for  each. 
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IV.  Analysis  of  Results 


The  multi-station  least  squares  data  fitting  procedure 
outlined  in  the  previous  section  has  been  applied  to  the 
simulated  data  scans  described  in  Section  II.  We  have  obtained 
estimates  of  the  position  and  velocity  uncertainties  as  a 
function  of  fence  crossing  longitude  and  altitude  for  each  of 
following  system  configurations: 

1)  .  Current  receiver  configuration  (six  sites)  with  phase 

and  Doppler  measurements  at  all  sites. 

2) .  Current  configuration  with  phase,  Doppler,  and  chirp 

measurements  at  all  sites. 

3)  .  Current  and  3  OOPS  locations  with  phase  and  Doppler 

measurements  at  all  sites. 

4)  .  Current  and  3  OOPS  locations  with  phase,  Doppler  and 

chirp  measurements  at  all  sites. 

5)  .  Current  and  3  OOPS  locations  with  phase  measurements 

at  current  sites  only,  Doppler  at  all  sites. 

6) .  Current  and  3  OOPS  locations  with  phase  measurements 

at  current  sites  only,  Doppler  and  chirp  at  all 

sites. 

For  each  of  these  configurations  we  have  obtained 
predictions  of  the  system  response  to  130  test  orbits  which 
sample  the  range  of  altitudes  and  fence  crossing  longitudes 
enumerated  in  Section  Ila.  The  assumed  orbital  inclination  was 
85°  and  the  assumed  radar  cross-section  was  1  meter2.  In 
general,  the  predicted  amplitudes  for  orbits  above  4000  km  were 
such  that  little  or  no  useful  data  was  generated,  so  we  exclude 
these  orbits  from  further  consideration. 

With  each  of  the  above  configurations  we  have  run  the  full 
grid  of  model  orbits  for  differe.it  Doppler  and,  where  included, 
chirp  accuracies.  The  range  of  Doppler  measurement  accuracy 
investigated  was  ±1  Hz  to  ±20  Hz,  while  the  chirp  accuracy  was 
varied  between  ±4  Hz/sec  and  ±20  Hz/sec. 

The  quantitative  results  from  our  model  runs  are  summarized 
in  Tables  la-12d.  These  Tables  list  for  each  case  the  expected 
total  (3-dimensional)  position  and  velocity  errors  as  a  function 
of  fence  crossing  position,  and  provide  a  breakdown  of  the 
velocity  errors  into  NAVSPASUR  great  circle  in-plane  (2-D)  and 
plane-normal  components.  It  is  clear  from  an  examination  of 


10 


these  Tables  that  the  measurement  accuracy  of  the  NAVSPASUR 
system  for  any  given  configuration  cannot  be  characterized  by  a 
single  number,  as  the  expected  error  depends  sensitively  on  the 
longitude  and  altitude  of  fence  crossing.  In  general,  all  cases 
studied  exhibit  their  best  performance  for  fence  crossing 
longitudes  between  85°W  and  100°W,  at  altitudes  between  200  and 
1600  km.  This  is  a  result  of  the  higher  density  of  receiver 
stations  to  the  east  of  the  Kickapoo  transmitter  site. 

We  provide  in  the  following  sections  a  detailed  analysis  of 
the  system  performance  for  the  various  options  examined. 

a).  Options  utilizing  the  current  receiver  configuration 

The  first  case  studied  corresponds  to  the  current  capability 
model  as  described  in  Table  5-3  of  the  System  Development  Plan, 
Phase  II  (hereafter  SDP)  .  For  this  case  we  assume  a  Doppler 
measurement  accuracy  of  *20  Hz  with  no  chirp.  In  discussions 
with  the  staff  at  NAVSPASUR  headquarters  we  have  been  advised 
that  the  phase-fitting  method  of  obtaining  the  Doppler  frequency 
yields  a  typical  accuracy  of  ±3-4  Hz.  However,  for  consistency 
with  the  SDP,  we  have  run  a  *20  Hz  model.  The  results  of  this 
model  are  summarized  in  Tables  la-d. 

An  examination  of  Table  la  indicates  that,  with  the 
exception  of  low  altitude  satellites  in  the  west,  the  present 
system  meets  the  position  determination  objective  of  *1.67  km 
over  CONUS  up  to  an  altitude  of  3000  km.  The  velocity 
performance  of  the  current  system,  however,  does  not  meet  the  *2 
km/min  requirement  at  any  fence  crossing  locations  with  the 
exception  of  low  altitude  satellites  at  90°  west  longitude  (see 
Table  lb).  A  comparison  of  in-plane  vs.  plane-normal  velocity 
performance,  listed  in  Tables  lc  and  Id  respectively, 
demonstrates  that  the  velocity  errors  are  dominated  by  the 
uncertainties  in  the  plane-normal  component.  An  analysis  of  the 
Doppler  measurement  equations  for  the  NAVSPASUR  system  shows 
that,  due  to  the  peculiar  geometry  of  the  current  system,  the 
received  Doppler  frequency  is  insensitive  to  the  plane-normal 
component  of  the  satellite  velocity,  so  that  improved  Doppler 
tracking  accuracy  from  in-plane  receivers  will  not  substantially 
affect  the  overall  velocity  accuracy.  We  include  a  full  analysis 
of  the  Doppler  measurement  equations  in  Appendix  A,  to  which  the 
interested  reader  is  referred. 

The  assesment  that  improved  Doppler  accuracy  will  not 
improve  overall  velocity  accuracy  is  borne  out  by  the  results  of 
our  second  model,  which  was  based  on  the  current  receiver 
configuration  with  an  assumed  a  Doppler  accuracy  of  *4  Hz.  While 
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the  improved  Doppler  accuracy  yields  in-plane  velocity 
measurements  which  are  a  factor  of  5  better  than  the  first  case 
(compare  Tables  lc  and  2c)  ,  there  is  no  net  improvement  in  the 
plane-normal  velocity  error  (Tables  Id  and  2d)  or  in  the 
resulting  total  velocity  error  (Tables  lb  and  2b) . 

Although  the  received  Doppler  frequency  is  insensitive  to 
plane-normal  velocity  changes,  the  chirp  does  contain  some 
information  on  plane-normal  velocity  (see  Appendix  A  for 
details) .  We  have  therefore  investigated  the  effect  of  utilizing 
chirp  information  in  the  fitting  process.  For  the  current 
receiver  configuration,  we  find  that  a  chirp  measurement  accuracy 
of  ±20  Hz/sec  combined  with  a  Doppler  measurement  accuracy  of  ±20 
Hz  yields  a  typical  velocity  improvement  of  2-4%  as  compared  to 
the  ±20  Hz  Doppler  measurement  alone  (ref.  Table  3b)  .  Clearly, 
if  chirp  is  to  yield  a  significant  improvement  to  the  velocity 
measurement  error,  chirp  measurement  accuracies  much  better  than 
±20  Hz/sec  are  required.  As  an  example,  if  the  chirp  measurement 
accuracy  is  improved  to  ±4  Hz/sec,  together  with  a  Doppler 
measurement  accuracy  of  ±1  Hz,  velocity  errors  are  reduced  by  30- 
50%  or  more  (see  Tables  4b, d) .  This  level  of  Doppler  and  chirp 
measurement  accuracy  represents  a  reasonable  limit  to  what  can  be 
achieved  with  the  current  receiver  system,  based  on  our  fits  to 
NAVSPASUR  phase4  data. 

It  is  apparent  from  Table  4b  that,  even  with  ±1  Hz  Doppler 
and  ±4  Hz/sec  chirp  capability,  the  current  system  cannot  achieve 
the  target  requirement  of  ±2  km/min  velocity  error  over  most  of 
the  fence.  This  limitation  results  from  the  planar  geometry  of 
the  current  system.  In  order  to  correct  this  deficiency,  one 
must  either  view  the  satellite  from  an  off-plane  location  while 
it  lies  in  the  NAVSPASUR  great  circle  plane,  or  illuminate  the 
ratellite  while  it  lies  outside  the  plane. 

b).  Configurations  including  out-of-plane  stations 

We  next  investigated  the  effect  of  adding  three  out-of-plane 
receiving  stations  to  the  current  system.  As  we  have  discussed 
in  Section  lie,  we  believe  this  to  be  the  minimum  number  of  OOPS 
required  to  give  adequate  coast-to-coas-t  performance.  Two 
different  types  of  OOPS  were  considered.  In  one  set  of  cases  we 
considered  the  effects  of  adding  OOPS  which  provided  both  phase 
and  Doppler  data,  while  in  the  other  case  we  examined  the  effect 
of  limiting  the  OOPS  to  Doppler  measurements  alone.  For  both 
cases  we  also  considered  the  impact  of  adding  chirp  measurement 
capability. 


i.  Doppler-only  OOPS 


In  general,  Doppler  measurements  made  from  receiving 
stations  well  removed  from  the  NAVSPASUR  great  circle  plane 
contain  significant  information  concerning  the  plane-normal 
velocity  (see  Appendix  A  for  details) .  For  systems  employing 
suitably  located  OOPS  receivers,  we  expect  the  Doppler-derived 
plane-normal  velocity  errors  to  be  comparable  to  the  in-plane 
errors,  and  we  further  expect  the  improvement  in  velocity 
determination  to  be  proportional  to  the  improvement  in  Doppler 
measurement  accuracy. 

Our  models  indicate  that  the  addition  of  three  Doppler-only 
OOPS  as  described  in  Section  lie  with  ±20  Hz  Doppler  accuracy 
would  yield  up  to  a  factor  of  three  improvement  in  velocity 
performance  over  the  current  configuration  with  the  same  Doppler 
accuracy  (ref.  Table  5b) .  As  expected,  this  improvement  comes 
primarily  from  an  enhanced  plane-normal  velocity  determination 
capability,  as  seen  from  a  comparison  of  Table  5c  with  Table  lc 
and  Table  5d  with  Table  Id.  While  the  improvement  is 
significant,  however,  the  velocity  errors  still  exceed  the  target 
requirement  over  much  of  the  fence. 

If  the  Doppler  accuracy  is  improved  to  ±10  Hz  the  velocity 
errors  are  reduced  by  an  additional  factor  of  ~2 ,  resulting  in 
errors  that  fall  within  the  target  requirement  over  most  of  the 
fence  at  altitudes  up  to  1300  km  (see  Table  6b) . 

In  contrast  with  the  results  obtained  for  the  system  based 
on  the  current  six  receivers,  the  addition  of  chirp  measurement 
capability  to  the  OOPS  configurations  yields  no  significant 
improvement  in  system  performance.  We  show  in  Tables  7a-d  the 
results  we  obtain  for  an  OOPS  configuration  with  ±10  Hz  Doppler 
accuracy  and  ±10  Hz/sec  chirp  accuracy.  These  results  are  almost 
identical  to  those  in  Tables  6a-d,  where  no  chirp  capability  was 
assumed. 

The  final  Doppler-only  OOPS  configuration  we  investigated 
assumed  a  Doppler  accuracy  of  ±1  Hz  and  a  chirp  accuracy  of  ±4 
Hz/sec  at  all  stations.  This  configuration  would  produce 
velocity  uncertainties  well  within  the  target  requirement  at 
virtually  all  fence  crossing  locations  up  to  an  altitude  of  4000 
km  (see  Tables  8a-d)  .  We  should  note  here  that  for  this 
configuration  the  positional  errors  do  exceed  the  target 
requirement  at  the  edges  of  CONUS  at  altitudes  above  3000  km. 
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ii.  Phase  and  Doppler  OOPS 


The  measurement  of  phase  offsets  between  spatially  separated 
antenna  arrays  at  the  OOPS  locations  can  provide  potentially 
dramatic  improvement  in  system  performance.  For  a  phase-capable 
OOPS  site  as  described  in  Section  lie  with  Doppler  accuracy  of 
*20  HZ,  we  obtain  velocity  errors  within  the  requirement  over 
most  of  CONUS  up  to  2500  km  altitude  (Table  9b) .  The  in-plane 
and  plane-normal  velocity  errors  for  this  configuration  are 
roughly  equal  (Tables  9c, d).  Moreover,  the  position  errors  are 
vastly  improved  over  the  systems  previously  discussed  (Table  9a) , 
with  errors  of  100  meters  or  less  over  most  of  the  fence  up  to 
4000  km  altitude. 

Improving  the  Doppler  accuracy  to  ±10  Hz  nets  an  additional 
velocity  improvement  of  5-15%  over  most  of  the  fence  (Tables  10a- 
d) ,  effectively  extending  the  altitude  limit  for  meeting  the 
velocity  requirement  to  3000  km.  Again,  we  find  that  adding 
chirp  capability  produces  no  significant  gains,  with  the 
exception  of  some  marginal  improvement  for  low  altitude  orbits 
near  the  edges  of  CONUS  (Tables  lla-d) . 

Finally,  we  analyzed  the  system  performance  for  an  assumed 
Doppler  accuracy  of  *1  Hz  and  a  chirp  accuracy  of  *4  Hz/sec.  The 
results,  shown  in  Tables  12a-d,  indicate  that  such  a  system  would 
be  capable  of  achieving  the  target  position  and  velocity 
requirements  up  to  an  altitude  of  4000  km  over  all  of  CONUS  with 
the  exception  of  100  km  orbits  at  the  extreme  edges  of  CONUS. 
Further,  this  system  could  provide  off-coast  (up  to  5°  longitude) 
performance  within  the  target  requirements  for  orbits  between  500 
and  4000  km.  Typical  velocity  errors  over  most  of  the  fence  are 
less  than  0.5  km/min,  while  typical  position  errors  are  less  than 
100  meters. 


V.  Identification  of  Further  Work 

Further  research  is  either  currently  underway  or  planned  in 
several  areas  related  to  the  work  described  in  this  report.  The 
highest  priority  research  effort  is  to  explore  the  possible  range 
of  number,  location,  and  operational  design  specifications  for 
the  OOPS.  The  current  work  demonstrates  that  the  required 
position  and  velocity  accuracies  can  be  obtained  with  3  OOPS  with 
the  performance  characteristics  assumed  in  this  report.  However, 
a  detailed  investigation  is  required  into  the  tradeoffs  between 
the  number  of  OOPS  sites  and  their  design  and  location  in  order 
to  optimize  system  performance  in  a  cost  effective  manner.  In 
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support  of  this  optimization  process  we  are  instituting  several 
enhancements  to  the  system  modelling  algorithms. 

As  we  have  have  shown  in  this  report,  the  uncertainties  in 
the  satellite  velocity  depend  sensitively  on  the  assumed  Doppler 
measurement  accuracy  at  the  OOPS  locations.  We  are  therefore 
analyzing  current  phase4  data  to  determine  the  level  of  Doppler 
and  chirp  accuracy  which  can  be  attained  with  the  present 
receiver  setup.  This  analysis  will  allow  us  to  generate  a  more 
realistic  Doppler/ chirp  error  model,  where  the  expected  errors  in 
these  parameters  depend  on  the  received  signal-to-noise  ratio  and 
the  number  of  useable  data  lines. 

We  are  also  working  to  develop  a  realistic  OOPS  gain  model, 
based  on  antenna  pattern  calculations  for  a  two-dimensional 
phased  array  of  multi-element  dipoles.  These  calculations  will 
allow  us  to  specify  the  OOPS  gain  as  a  function  of  azimuth  and 
elevation  to  the  target  satellite  as  opposed  to  a  constant-gain 
model  as  assumed  for  this  report. 

When  the  above  enhancements  are  incorporated  into  the  sytem 
model,  we  will  proceed  with  the  determination  of  the  optimum 
system  configuration  and  the  development  of  design  specifications 
for  the  OOPS  receivers. 

In  addition  to  the  above-mentioned  efforts,  we  recognize  the 
need  to  extend  this  work  to  include  performance  estimates  for 
highly  elliptic  satellite  orbits,  and  to  investigate  the 
performance  of  both  the  half-Doppler  and  quarter-Doppler  modes. 

V.  Summary 

We  have  developed  a  set  of  algorithms  which  allow  us  to 
accurately  model  the  full -Doppler  mode  performance  of  the  current 
NAVSPASUR  system,  as  well  as  the  proposed  enhancements  outlined 
in  the  NAVSPASUR  System  Development  Plan,  Phase  II.  The 
enhancements  which  have  been  studied  to  date  fall  into  two 
categories:  1)  systems  based  on  the  current  six  receiver  sites 
but  requiring  additional  and/or  improved  signal  processing,  and 
2)  systems  requiring  the  addition  of  three  new  out-of-plane 
receiving  stations,  either  with  or  without  enhanced  signal 
processing  capability. 

The  results  of  our  analysis  indicate  that  none  of  the 
enhancements  which  are  based  solely  on  the  current  six  receiving 
stations  are  capable  of  meeting  the  stated  velocity  accuracy  goal 
of  *2  km/min  over  all  of  CONUS.  With  the  addition  of  three  OOPS 
with  both  phase  and  Doppler/ chirp  measurement  capability,  it  is 


possible  to  meet  the  velocity  accuracy  requirement  over  CONUS  up 
to  an  altitude  of  4000  km.  If  the  three  OOPS  are  configured  for 
Doppler/ chirp  measurement  only,  velocity  errors  within  the  stated 
specification  may  still  be  possible  over  CONUS  up  to  4000  km  if 
Doppler/ chirp  measurement  accuracies  of  *lHz  and  *4Hz/sec, 
respectively,  can  be  achieved. 

We  have  demonstrated  in  this  work  that  there  exists  at  least 
one  system  enhancement  geometry  which  is  capable  of  providing 
single-pass  position  and  velocity  accuracies  meeting  the  target 
requirements  specified  in  the  NAVSPASUR  System  Development  Plan, 
Phase  I.  The  results  presented  herein  indicate  that  these 
requirements  can  be  met  with  an  expanded  system  which  would 
include  three  OOPS  receiver  sites  utilizing  the  current 
transmitters  for  satellite  illumination.  While  our  results  show 
that  the  accuracy  requirements  can  be  met  with  the  system 
enhancements  proposed  herein,  additional  work  must  be  done  to 
fully  explore  the  range  of  possible  design  and  location  options 
for  the  OOPS.  We  are  currently  investigating  the  technical  issues 
relating  to  the  OOPS  design  and  location,  and  will  present  our 
findings  in  an  upcoming  report. 
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7IGUR2  1:  Map  depicting  the  geographic  locations  of  the  current 
NAVSPASUR  transmitter  and  receiver  stations,  as  well  as  the 
nominal  positions  of  the  three  OOPS  receiver  sites  used  in  this 
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APPENDIX  A 


This  appendix  presents  a  brief  discussion  of  how  Doppler  and 
Doppler  rate  (chirp)  measurements  can  be  used  to  improve  the 
velocity  and  position  determination  in  the  NAVSPASUR  system. 

For  a  bi-static  radar,  the  observed  Doppler  shift  frequency 
vD  is  given  by 


vD(t) 


(A.  1) 


where  v0  is  the  NAVSPASUR  transmit  frequency  (216.98  MHz),  c  is 
the  velocity  of  light,  vs  is  the  velocity  vector  of  the  target 
satellite,  and  rls,  rrs  are  the  unit  vectors  to  the  satellite  from 
the  transmitter  and  receiver,  respectively.  Since  for  all  orbits 
of  interest  (vs  ♦  rts)  «  c,  equation  (A.l)  reduces  to 

M*)  *  ~  ^  {  vs  •  [rts  +  rrs]  ).  (A.  2) 


Adopting  an  earth-centered  cartesian  coordinate  system  where 
x,  y  are  in  the  NAVSPASUR  great  circle  plane  and  z  is  normal  to 
the  plane,  equation  (A. 2)  becomes 


VD  (^0) 


V0  f  vx(xs-x,)  +  vy(ys-yi)  +  vx(zs-zt) 

c  1  ej 


vx(xs-xr)  +  vy(ys-yr)  +  vx(zs-zr) 


(A. 3) 


where 


lrtsl  =  [  (xs-xt)2  +  (ys-yt)2  +  (zs-zt)2  \W2 

1/2 

lrrsl  =  j  (xs-xr)2  +  (ys-yr)2  +  (zs-zr)2  ]  . 
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Here  vx,  vy,  vz  are  the  components  of  the  satellite  velocity  at 
time  t0  and  xs,  ys,  zs,  are  the  satellite  position  coordinates  at 
t0.  The  transmitter  and  receiver  coordinates  are  given  by  xt,  yt, 
zt  and  xr,  yf,  zr  respectively.  For  NAVSPASUR  observations  t0, 
the  time  at  which  the  satellite  position  and  velocity  are 
determined,  is  nominally  the  time  of  fence  crossing. 


Equation  (A. 3)  provides  a  means  to  examine  the  relative 
information  content  of  the  Doppler  frequency  in  light  of  the 
special  geometry  of  the  current  NAVSPASUR  system.  Specifically, 
the  sensitivity  of  the  Doppler  frequency  to  any  satellite 
position  or  velocity  component  is  given  by  the  partial  derivative 
of  equation  (A. 3),  taken  with  respect  to  that  component.  Taking 
the  partial  derivatives  with  respect  to  the  velocity  components, 
for  example,  yields  the  following  equations: 


aVD 

_  _  uo 

(xs-xt)  (  (xs-xr) 

3VX 

>:>0  C 

®„i  »„i  1 

3wo 

[<V*.>  ,  <Vy,>  1 

3vy 

t-t  c 

IM  1 

3vd 

[<*,•*,)  ,  (Zs"Zt>  1 

3Vz 

t=lo  C 

pj  p„l 

In  deriving 

the  above  expressions  we 

partial  derivatives  at  time 

tQ  in  order  to  si 

(A. 4) 


(A. 5) 


(A. 6) 


expressions  and  to  clarify  the  functional  dependences,  over  the 
short  period  during  which  the  satellite  is  illuminated,  the  time 
dependence  of  the  partial  derivatives  is  weak,  so  that  our 
analysis  is  not  seriously  limited  by  this  simplification. 


In  the  coordinate  system  we  have  selected,  both  (zs-zt)  and 
(zs-zr)  are  much  smaller  than  |?  |  and  |?  |  respectively,  since  all 
of  the  receiver  and  transmitter  sites  lie  within  the  NAVSPASUR 
plane,  and  the  satellite  is  only  illuminated  within  a  small 
distance  of  the  plane.  This  implies  that  the  Doppler  frequency 
is  relatively  insensitive  to  vz  (t.e.,  3vD/3vz  is  small  compared  to 
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the  remaining  velocity  partial  derivatives) .  Note  that  the 
corresponding  case  is  almost  never  true  for  the  two  in-plane 
velocities  because  of  the  bi-static  nature  of  the  NAVSPASUR 
system.  In  other  words,  if,  for  example,  (xs-xt)  ~  (zs-z(),  then 
(xs-xr)  »  (zs-zr),  since  the  transmitter  and  receiver  sites  are 
well  separated.  For  this  reason  the  Doppler  frequency  will 
always  contain  significant  information  on  both  in-plane 
velocities.  For  typical  distances  encountered  with  NAVSPASUR 
observations,  the  two  in-plane  partial  derivatives  are  of  order 
10  Hz  per  km/min,  so  that  a  Doppler  measurement  accuracy  of  ±10 
Hz  or  better  can  constrain  the  in-plane  velocity  within  the  FPOD 
requirements.  Typical  values  of  the  plane-normal  partial 
derivatives,  however,  are  of  order  0.01  Hz  per  km/min,  so  that 
the  plane-normal  velocity  cannot  be  adequately  constrained  for 
any  reasonably  achievable  Doppler  accuracy. 


One  can  see  from  equation  (A. 6)  that  Doppler  measurements 
from  an  OOPS  receiver  would  provide  enhanced  sensitivity  to  the 
plane-normal  velocity,  since  for  stations  well  removed  from  the 
great  circle  plane  (zs-zr)  would  be  comparable  to  |rrs|. 


In  addition  to  containing  information  on  the  in-plane 
velocity  components,  the  Doppler  frequency  also  contains  some 
information  on  all  three  position  components.  Taking  the  partial 
derivatives  of  equation  (A. 3)  with  respect  to  the  position 
components  and  evaluating  at  time  t0,  we  have 


3vd  I 

1 

„  V0  V*  I 

f  l 

4- 

l 

3XS  1 

uto 

c  1 

■ 

i 

3V0  I 

_  v°  \  1 

l 

+ 

1 

sy  1 

l=‘o 

C 

i  |f,  | 

|f  1 

S 

1  ts' 

1  rs1 

3vd  I 

_  VQ  Vz  I 

■  1 

+ 

1 

3ZS  1 

uto 

c  : 

i 

(A. 7) 


(A. 8) 


(A. 9) 


The  relative  sensitivity  of  the  Doppler  to  each  of  the  position 
coordinates  depends  on  the  specifics  of  the  satellite  orbit,  but, 
in  general,  the  three  velocity  components  are  comparable,  so  that 
the  sensitivity  to  each  position  coordinate  is  comparable. 
Typical  values  for  the  position  partial  derivatives  are  of  order 
1  Hz  per  km,  so  that  Doppler  information  is  not  expected  to 
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significantly  improve  upon  the  position  determinations  derived 
from  the  phase  data. 

In  the  case  of  chirp  measurements,  ve  have,  taking  the  time 
derivative  of  equation  (A. 2), 

*o  "  -  J  (  **  *  +  +  **  '  dE  }  (A-10) 

where  the  satellite  acceleration,  as,  is  given  by 

-  GM*  . 

as  -  -  R.  (A.  11) 

Here,  R  is  the  satellite  position  vector  (earth-centered) ,  G  is 
the  gravitational  constant,  and  Mg  is  the  mass  of  the  earth.  For 
typical  values  of  the  distances  encountered  with  NAVSPASUR 
satellite  observations,  the  first  term  of  equation  (A. 10)  is 


GMg  V„  . 

gjF  R  •  [ru  +  rrs]  ^  <  7  HZ  sec-1.  (A.  12) 

Since  the  satellite  can  only  be  detected  if  it  is  above  the 
horizon  of  both  the  transmitter  and  receiver,  both  dot  products 
in  equation  (A. 12)  are  positive,  so  that  the  contribution  of  the 
acceleration  term  to  the  observed  chirp  is  positive,  and,  as  we 
show  below,  is  a  small  part  of  the  total  chirp.  Further,  the 
contribution  of  this  term  depends  only  on  the  satellite  position 
and  may  be  calculated  from  the  position,  which  is  normally  well 
determined  by  the  data. 

The  second  term  of  equation  (A. 10)  reduces  to 


(A. 13) 
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The  terms  in  the  square  brackets  in  equation  (A. 13)  are 
always  greater  than  or  equal  to  zero,  since  |$s|  is  greater  than  or 
equal  to  both  |vs  •  rt$|  and  |tfs  •  rrs|.  Further,  due  to  the  bi-static 
nature  of  the  system,  if  one  of  the  terms  in  the  square  brackets 
is  zero,  the  other  is  not.  Therefore,  the  net  chirp  due  to  the 
second  term  in  equation  (A. 10)  is  negative  for  all  satellite 
passes.  For  typical  NAVSPASUR  satellite  observations,  the  chirp 
contribution  of  the  second  term  is  between  -20  and  -100  Hz/sec. 

As  in  the  case  of  the  Doppler  frequency,  we  can  examine  the 
information  content  of  the  chirp  by  taking  the  partial 
derivatives  of  equation  (A. 10)  with  respect  to  the  position  and 
velocity  components.  Taking  the  velocity  partials,  again  at  time 
t0,  we  have 


*1  I 

3v*  H 


(A. 14) 


(A. 15) 


(A. 16) 


Equations  (A. 14)  -  (A. 16)  demonstrate  that  the  chirp  will  in 
general  contain  information  on  all  three  velocity  components. 
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Because  both  (zs-zt)  «  !rtsl  and  (zs-zr)  «  lrrsl,  equation  (A.  16) 
is  approximately 


3vn  2  V,  Vg 

- - 


1  +  1 


c  1 


(A. 17) 


The  sensitivity  of  the  chirp  to  the  plane-normal  velocity  is 
thus  proportional  to  the  plane-normal  velocity.  In  the  case  of 
the  in-plane  components,  the  functional  dependence  of  the  partial 
derivatives  is  somewhat  more  complicated.  However,  one  can 
deduce  from  an  examination  of  equations  (A. 14)  and  (A. 15)  that 
the  contribution  of  the  in-plane  terms  is  generally  of  the  same 
order  of  magnitude  as  for  the  plane-normal  term.  This  result  is 
in  contrast  to  the  result  obtained  for  the  Doppler  measurement, 
where  the  plane-normal  velocity  sensitivity  is  negligible.  For 
typical  NAVSPASUR  observations,  the  partial  derivatives  of  the 
chirp  with  respect  to  the  velocity  components  are  of  order  0.05 
Hz/second  per  km/min.  In  general,  then,  chirp  measurements  do 
not  provide  as  stringent  a  constraint  on  the  velocity  as  the 
Doppler  measurements. 


Summary 


We  ha /e  shown  that,  in  general,  both  the  Doppler  and  chirp 
measurements  contain  information  on  all  satellite  position  and 
velocity  components.  Doppler  measurements  from  in-plane 
receiving  sites  are  far  more  sensitive  to  in-plane  velocity  than 
to  plane-normal  velocity.  Given  the  typical  magnitude  of  the 
Doppler  partial  derivatives  with  respect  to  velocity,  we  expect 
that  Doppler  measurement  accuracies  of  ±10  Hz  can  provide  in¬ 
plane  velocity  determinations  of  order  ±1  km/min,  much  better 
than  that  attainable  from  direction  cosine  rates  derived  from  the 
phase  data.  However,  for  the  current  system,  the  Doppler 
frequency  cannot  constrain  plane-normal  velocity  nearly  as  well 
as  the  ph7ise  data.  In  order  to  achieve  plane-normal  velocity 
accuracies  comparable  to  the  in-plane  velocity  accuracy,  either 
the  target  satellite  must  be  illuminated  well  off  the  NAVSPASUR 
great  circle  plane,  or  out-of-plane  receiver  stations  must  be 
added ,  or  both . 

In  the  case  of  position  information,  Doppler  measurements 
cannot  substantially  improve  upon  the  accuracy  attainable  from 
the  direction  cosines  derived  from  the  phase  data. 
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Finally,  while  chirp  measurements  yield  in-plane  and  plane- 
normal  velocity  accuracies  which  are  comparable,  in  neither  case 
is  the  sensitivity  sufficient  to  yield  significant  improvement 
over  the  accuracies  available  from  either  the  phase  or  Doppler 
information. 
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longituda.  Tha  noainal  latituda  is  obtainad  (roa  a  cubic  aquation 
fit  to  tha  positions  of  tha  racaiving  stations.  This  aathod 
pcovidas  sufficiant  accuracy  to  proparly  position  tha  sataliita 
within  a  faw  kiloaatara  of  tha  fanca. 


Convert  tho  satellite  position  and  velocity  to  NAVSPASUB  great 
circle  coordinates. 
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Wt  will  obtain  tha  gaocantric  coordinates  o f  tha  satallita  froa  tha 
gaodatic  latitude,  longitude  and  altitude  via  FUND.  OF  ASTR.  page  9$ 
aquations  (2.8-7)  and  (2.8-8).  Tha  first  tarns  in  aquations  (2.8-7) 
depend  only  on  latitude,  so  calculate  than  now. 
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NONROTATING  earth. 
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rind  a  aaxiaua  of  55  data  linas  with  aaplitudas  abova  tha  cutoff. 
Bagin  tha  calculations  TPRED  -  40*TINCR  saconds .  If  55  linas  ata 
not  found  aftar  110  trias  (2  saconds),  giva  up. 
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Calculate  the  power  at  the  satellite  based  on  the  transmitter  power, 
antenna  gain  in  the  direction  of  the  satellite,  and  t r ansmi t te r- to- 
satellite  distance. 
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rWHP .  Tha  correction  has  baan  zodallad  as  an  additiva  ten  which  is 
proportional  to  tha  antanna  array  langth  and  invarsaly  proportional 
to  tha  TX- s atallita  distanca.  Tha  factor  of  3  was  datarainad 
aapiricaXIy  by  coaparison  of  raal  and  siaulatad  data. 
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21.  Calculate  the  Doppler  residual  (observed  -  calculated). 
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This  is  th«  loopback  point  for  iterating  tha  solution. 


through  tha  scans  and  th#  antannas  within  a  scan, 
ork  is  all  dona  in  tha  subroutinas.  Tha  subroutinas  MCALC 
:0CRIV  ara  problaa  spacific  i.a.,  tha  natura  of  tha  obsarvablas 
ba  considarad.  Tha  raaaining  subroutinas  ara  ganaral.  Thay 
sad  to  sat-up  and  invart  tha  aatrieas  and  display  tha  rasults. 
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Calculate  the  difference  in  path  lengths  fro»  the  satellite  to  the  JNth 
antenna  and  the  reference  antenna,  respectively.  Express  the  result  in 
wavelengths  at  the  Doppler-shifted  receive  frequency. 
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